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Deep Chandra observations (53 ks, ACIS-S3) of NGC 3077, a starburst dwarf 
galaxy in the M81 triplet, resolve the X-ray emission from several supershells. 
The emission is brightest in the cavities defined by expanding shells detected 
previously in Ha emission (Martin 1998). Thermal emission models fitted to the 
data imply temperatures ranging from ~ 1.3 to 4.9 x 10 6 K and indicate that the 
strongest absorption is coincident with the densest clouds traced by CO emission. 
The fitted emission measures give pressures of P/k « 10 5 ~ 6 £~°- 5 /~ ' 5 K cm~ 3 
(£: metallicity of the hot gas in solar units, /„: volume filling factor). Despite 
these high pressures, the radial density profile of the hot gas is not as steep 
as that expected in a freely expanding wind (e.g., as seen in the neighboring 
starburst galaxy M82) implying that the hot gas is still confined by the Ha 
shells. The chaotic dynamical state of NGC 3077 undermines reliable estimates 
of the escape velocity. The more relevant quantity for the ultimate fate of the 
outflow is probably the gas density in the rich intragroup medium. Based on the 
H I distribution of NGC 3077 and a connected tidal tail we argue that the wind 
has the potential to leave the gravitational well of NGC 3077 to the north but not 
to the south. The total 0.3 — 6.0 keV X-ray luminosity is ~ 2 — 5 x 10 39 erg s _1 
(depending on the selected thermal plasma model). Most (~ 85%) of the X-ray 
luminosity in NGC 3077 comes from the hot interstellar gas; the remainder comes 
from six X-ray point sources. In spite of previous claims to the contrary, we do 
not find X-ray emission originating from the prominent tidal tail near NGC 3077. 

Subject headings: galaxies: individual (NGC 3077) — galaxies: dwarf — galaxies: 
ISM — galaxies: irregular — galaxies: starburst — galaxies: evolution — ISM: 
bubbles — ISM: structure — X-rays 

1. Introduction 

In lower mass galaxies, violent star formation activity is predicted to drive interstellar 
gas into either the dwarf galaxy's halo or the intergalactic medium (IGM; Dekel & Silk 1986; 
Mac Low & Ferrara 1999; Silich & Tenorio-Tagle 2001). Cold dark matter cosmologies predict 
that dwarf galaxies began forming before more massive galaxies could be assembled. The fate 
of the dwarf galaxy winds may therefore have influenced the formation of the larger galaxies. 
It is not well understood, however, whether the properties of the surrounding intragroup gas 
or the gravitational potential dictate the fate of such galactic winds. Although intergalactic 
gas densities are much lower at the present epoch, detailed examination of the wind/group 
interaction in the local universe should improve our understanding of it at all cosmic epochs. 
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It has long been known through deep observations in the Ha line that galactic outflows 
exist, but relatively little is known about the hot phase of the gas (T ~ 10 6 ~ 7 K) which ac- 
tually drives the outflow. Early studies with EINSTEIN and ROSAT of extremely luminous 
starburst galaxies revealed the presence of hot winds, sometimes reaching distances above 
the plane of several kpc (e.g., M82: Strickland, Ponman, & Stevens 1997 in the following 
S97; Lehnert, Heckman, & Weaver 1999, NGC253: Strickland, Heckman, Weaver, & Dahlem 
2000; Strickland et al. 2002; Pietsch et al. 2001). Measurements of the temperature and den- 
sity in the hot winds were critical for estimating their power requirements and establishing 
the mechanical power from supernovae in the starburst as the energy source (see, e.g, the 
reviews of Heckman, Lehnert, & Armus 1993; Heckman 2002). The latest generation of 
X-ray observatories, Chandra and XMM-Newton, allow similar measurements of smaller, 
fainter galaxies due to their larger collecting areas. The angular resolution of Chandra, ~ 1", 
also allows a) the direct identification of discrete X-ray sources and b) spectral imaging of 
extended gas structures on spatial scales attained in optical and radio maps. Such detailed 
information turns out to be crucial for establishing a physical model of the rather complex 
integrated X-ray spectrum. 

In this paper we present Chandra X-ray observations of the dwarf starburst galaxy 
NGC 3077. We chose this particular dwarf galaxy for our study because NGC 3077 shows 
prominent expanding Ha supershells (Martin 1998, in the following M98) and is a member 
of the gas-rich M81 group. Interaction with M81 and M82 is believed to have triggered star- 
burst activity in both M 82 and NGC 3077. Tidal tails in the M 81 triplet allow to investigate 
environmental effects on the development of superwinds. Various studies of NGC 3077 have 
been published in the past (e.g., optical studies: Barbieri, Bertola, & di Tullio 1974, atomic 
and molecular gas: Walter, Weiss, Martin, & Scoville 2002, in the following W02; Meier, 
Turner, & Beck 2001). Previous ROSAT X-ray observations indicating the presence of dif- 
fuse X-ray emission in NGC 3077 were presented by Bi, Arp, & Zimmermann (1994). The 
interpretation of these observations, however, are complicated due to point source confusion. 

We adopt a distance to NGC 3077 of 3.6 Mpc (as obtained for M81, Freedman et al. 
1994), consistent with more recent studies (cf. Sakai & Madore 2001; Karachentsev et al. 
2002). In Sect. 2 we describe the Chandra observations, the data reduction, and the results. 
A discussion is presented in Sect. 3 followed by conclusions and a summary in Sect. 4. 

2. Observations, Data Analysis, and Results 

We observed NGC 3077 with the Chandra X-ray Observatory on 2001 March 07 and 
08 for 53.4 ks with the back-illuminated Advanced CCD Imaging Spectrometer (ACIS) S3 
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CCD chip (pixel size: 0'/49 x 0'.'49, energy resolution ~ 120 eV at IkeV; observation ID: 
2076, Seq.No. 600210). ACIS was at a temperature of — 120° C during these observations. 
The satellite telemetry was processed at the Chandra X-ray Center (CXC) with the Stan- 
dard Data Processing (SDP) system to correct for the motion of the satellite and to apply 
cosmic ray rejection. Since the data were processed before 2001 September 13, we applied 
a new gain map to accommodate the latest Order-sorting/integrated Probability Tables 
(OSIP). The data products were then analyzed with the CXC Chandra Interactive Analysis 
of Observations (CIAO) software version 2.2. CIAO Data Model tools were used for data 
manipulation, and Sherpa v2.2 for the fitting of the ACIS spectra. CALDB v2.9 of the ACIS 
CCD calibration files were applied throughout the analysis. We constructed a lightcurve by 
binning the dataset to a time resolution of 200 s. No significant flaring of the background 
was observed. Bad pixels of all images were screened out by applying the bad pixel mask 
provided by the CXC. We checked the Chandra pointing during our observations by verifying 
the position of the bright double star HD 86677. Its position falls within 1" of the measured 
coordinates listed in the HIPPARCOS Catalog (Perryman et al. 1997). 



2.1. Global Properties 

We defined four broad energy bands for the analysis of the X-ray emission of NGC 3077 
in an attempt to separate strong oxygen lines (soft band: 0.3 keV^ E ^0.7 keV) from major 
iron-L line complexes (medium band: 0.7keV^ E ^ 1.1 keV) of a hot thermal plasma. Line 
emission does not play an important role for higher energies. The definition of the third 
band (hard band: 1.1 keV^ E ^ 6.0 keV) therefore provides a good representation of the 
continuum emission. As photoelectric absorption (cross section a) is a function of energy 
E (a oc E~ 3 ), it is strongest in the soft band. The hard band, in contrast, is a better 
tracer for temperature. We also defined a total band (0.3 keV^ E ^ 6.0 keV) to discuss 
the global characteristics of NGC 3077. The resulting images were adaptively smoothed 
with the CIAO task csmooth, using a fast Fourier algorithm and a minimum and maximum 
significance S/N level of 2 and 3, respectively. The levels were chosen this way in order to 
give the best representation of features which are visible in the raw data set. Adaptively 
smoothed images, however, were not used for any quantitative analysis. 

Exposure maps were created for the soft, the medium, and the hard band based on their 
central energies. Since the intensity of the diffuse emission peaks at about 1 keV, we decided 
to use the exposure map of the medium band for the total band, too. The smoothing 
kernel of the adaptively smoothed images changes size for each pixel, and exposure map 
correction must account for this effect. To do so, we smoothed the exposure maps by using 
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the smoothing scales which were previously computed for the adaptively smoothed images. 
All images were finally corrected by dividing by the appropriate exposure map. 

The distribution of the X-ray emission in the different bands are shown in Fig. 1 - 
Fig. 2 a is a three color composite of the same images. Point sources contribute ~27% to 
the total count rate of NGC3077 (see Sect. 2.2.1) and are labeled in Fig. 1; the remaining 
63% is diffuse emission. The diffuse X-ray emission of NGC3077 has a size of about 1' 
(corresponding to 1 kpc). The soft emission is slightly more extended than the medium and 
hard emission (see discussion in Sect. 2.3), especially to the east of NGC3077 (see Fig. 1). 

We constructed two hardness ratios from the bands defined above: 



HR1 is more sensitive to photoelectric absorption than HR2, while HR2 can be used to 
trace temperature variations. The advantage of defining the hardness ratios via three bands 
compared to a more classical two band definition is that the values are still defined for all 
sources even if virtually all emission emerges from within one band. 

In Fig. 3 we present the hardness ratio images of NGC 3077 based on the adaptively 
smoothed broad band images. Due to the non-linearity of the csmooth algorithm, deviations 
from the real distributions may appear. Still, this representation gives a good overview on 
the global characteristics of the hardness ratios within NGC 3077. We find that the point 
source S4 is the softest source and S2 and S3 are the hardest. The hardest diffuse X-ray 
emission is found near the center of NGC 3077, close to SI, S2, and S3. To the north and 
west of NGC 3077, the X-ray emission is brighter in the medium band as compared to the 
soft band. 

Photoelectric absorption of X-ray emission leads to a hardening of the spectrum as it 
is most effective at low photon energies (oc E~ 3 ). One would expect that this effect happens 
where cold gas is located in front of the X-ray emission. To test this, we compare the HR1 
image to the H I and the CO distributions of NGC 3077 in Fig. 4 (data taken from W02; see 
also Fig. 2c). The dense molecular clouds are located where the center of NGC 3077 exhibits 
the hardest diffuse X-ray emission. Moreover, the H I emission follows nicely the HR1=— 0.3 
contour. In summary, the X-ray absorption features seen in NGC 3077 can potentially be 
explained by the presence of atomic and molecular gas. 



HR1 



HR2 



(Soft - Medium - Hard) /(Soft + Medium + Hard) 
(Soft + Medium - Hard) /(Soft + Medium + Hard) 



(1) 
(2) 
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2.2. Discrete Sources 

2.2.1. Source Detection 



In order to separate point sources from diffuse emission, we applied the wavelet source 
detection algorithm wavdetect to the spatially unbinned, total band image. The point spread 
function (PSF) degrades at large off-axis angles and we accounted for this effect by using 
0"5, l'.'O, and 2"0 wavelet radii. Using this method, 42 sources were detected on the ACIS- 
S3 chip; four of the detections are located within the diffuse X-ray and Ha emission of 
NGC 3077. To cross-check this result, we ran the alternative sliding cell celldetect algorithm 
on the same dataset. Regarding the varying shape of the PSF, the signal to noise (S/N) 
thresholds were chosen to be 1.9, 2.1, and 3.2, where the first threshold was applied within the 
inner 1' radius centered on the aimpoint. The second and third thresholds were used within 
radii r of 1' < r < 2'5 and 2'.5 < r < 3'5. Beyond 3f5, the PSF degrades substantially and 
detections were generally disregarded (for a justification see Martin, Kobulnicky, & Heckman 
2002). This led to a total detection of 28 sources, where again four belong to NGC 3077. 
No additional sources were detected by celldetect. Sources with low S/N were more likely 
to be detected by wavdetect than by celldetect. In Table 1 we list the detections within the 
ACIS-S3 field of view (excluding those located within the optical extent of NGC 3077; they 
are listed separately in Table 2). The source numbers (column 1) are the same as in Fig. 5. 
Column (2) lists the detection algorithm that detected the source (w: wavdetect, c: celldetect, 
x: the distance to the aimpoint is > 3f5; celldetect was not applied). The net source counts 
in the different bands are given in columns (5) to (8) and the hardness ratios HR1 and HR2 
in columns (9) and (10). A star at the source number denotes an optical counterpart from 
an R c -band image within a search radius of 10" (R c and B-band images obtained by us at 
the Calar Alto 2.2m telescope with integration times of 600 s and 900 s, respectively). This 
radius was chosen to account for seeing and positional uncertainties of the ground based 
data as well as for the extents of the tentative optical counterparts. 



2.2.2. Point sources in NGC 3077 



As mentioned above, both source detection algorithms detected four discrete sources 
within the diffuse X-ray and Ha emission of NGC 3077 which supposedly belong to this 
galaxy. From a visual inspection of the total band image, we decided to add a source close to 
the center of NGC 3077 to our list (S5, see Fig. 1). One detection exhibits two intensity peaks 
separated by 1" which is about the size of the ACIS-S3 PSF at IkeV (FWHM[PSF] « 1"1). 
We decided to split this detection into two sources (S2 and S3) for separate analysis. To check 
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whether or not the discrete sources are indeed point-like, we constructed a more accurate 
PSF for each source position, given the intensity-weighted mean energy of the corresponding 
source spectrum. Except for S5 which appears to be slightly extended (deconvolved size: 
l'/l x l'/3), all other sources are indeed point sources. 

Spectra of the point sources were extracted from apertures which were chosen to be three 
times larger than the corresponding PSF, except for S2 and S3; in these cases the spectra 
were extracted from rectangular boxes 3'/0 x 2'.'3 in size in an attempt to separate both 
objects (see Fig. 1). The corresponding Redistribution Matrix Files (RMFs) and Auxiliary 
Response Files (ARFs) were subsequently created for each source position on the CCD to 
calibrate the corresponding spectra. The spectral channels were not binned for any spectral 
fitting analysis in order to preserve the maximum of statistical information. One should keep 
in mind, however, that the spectral resolution of ACIS-S3 is about 120 eV at IkeV which 
corresponds to about eight pulse invariant (PI) energy channels. 

Each spectrum is composed of the real source spectrum and an energy-dependent offset 
which comprises instrumental and extragalactic backgrounds as well as diffuse X-ray emis- 
sion from NGC3077 itself. We checked the intensity of these offsets by extracting source 
counts in regions with an X-ray contribution similar to that expected at the source positions. 
This was obtained within apertures of the same size as those used for the extraction of the 
point source spectra. As only 1-3 counts of background were detected within each aperture, 
we decided to neglect the background contamination of the X-ray spectra of the discrete 
point source population. 

The final point source spectra can be divided into three types: a) flat spectrum sources 
with counts all across the ACIS energy range (sources S2 and S3), b) source spectra which 
peak at energies ~0.8-1.2keV (SI, S5, and S6), and c) a spectrum (S4), which peaks at 
~0.6keV and virtually shows no emission above 0.8 keV. 

The spectra were modeled within the energy range of 0.3 keV^ E ^ 8.0 keV (neglecting 
data with bad quality flags). We fit three different models to each spectrum: a Raymond- 
Smith collisional thermal plasma (RS, Raymond & Smith 1977, and updates 4 ), a power law 
(PL) and a black body (BB) model. All source models were coupled with the photoelec- 
tric absorption model from Balucinska-Church & McCammon (1992). Following Heckman 
(1980) and Martin (1997), we assumed solar metallicity for NGC3077. This value is in good 
agreement with the global optical spectrum of NGC 3077 (Kennicutt 1992) as applied to the 
oxygen abundance estimation method recently developed by Kewley & Dopita (2002). The 
solar element mixture, which is implemented in the Sherpa code was adopted from Anders 



see the XSPEC manual: http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/manual.html 
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& Grevesse (1989). The metallicity of NGC3077 is similar to that of our Galaxy. For this 
reason, it is sufficient to model the photoelectric absorption (due to the Galactic interstellar 
medium [ISM] as well as to material within NGC3077) by a single component, with the 
restriction that the derived column densities must exceed the Galactic value. From the Lei- 
den/Dwingeloo Survey (Hartman & Burton 1997) we measure a Galactic H I column density 
of 4 x 10 20 cm- 2 at the position of NGC3077. 

We used the Cash maximum likelihood statistics for best fit estimation, which is more 
accurate than the \ 2 statistics in the case of low number statistics. The optimization method 
was a Monte Carlo algorithm with 64 random starting values within a physically sensible 
parameter space. Note that a lower number of starting values provided unstable results. For 
each of these initial parameters, the single-shot Powell optimization method locally max- 
imized the likelihood (Monte-Powell within Sherpa). The fits with the supposedly global 
maximum likelihood were finally confirmed by a Levenberg-Marquardt optimization method 
algorithm. The best fit parameters for each source and model are shown in Table 2, where 
Nh is the absorbing column density, T the temperature, 7 the power law photon index, 
F x bs the absorbed flux, Fx the unabsorbed flux, and Lx the unabsorbed X-ray luminosity. 
Unfortunately, the surfaces of constant likelihood in the parameter spaces especially of ab- 
sorbed thermal plasma models are not well-behaved, i.e., they are not multi-dimensional 
paraboloids. The use of any standard error estimation algorithm like Sherpa's uncertainty, 
or covariance will therefore inevitably fail. For this reason, errors of the individual fits were 
obtained by using the "frequent ist's method" (see, e.g., Freeman, Doe, & Siemiginowska 
2001). For each set of best fit parameters we simulated 500 spectra with the Poisson noise 
distribution of the observations. Each of these simulations were re-fitted using the best fit 
values as initial values for a Powell optimization method. The resulting values for all pa- 
rameters were used to construct a sampling statistic and the inner 68% of the distribution 
obtained are considered to be good estimates for the 1 a errors. Errors for fluxes and lumi- 
nosities were calculated on the basis of extreme parameter values. The most plausible best 
fits are displayed in Figs. 6 and 7 (see Sect. 3.1 for a discussion) as well as the corresponding 
confidence regions in the iV^-T/7 planes. Unfortunately, the limited S/N does not allow to 
give preference to any specific model based on a determination of a goodness-of-fit indicator. 

We checked all discrete point sources for signs of variability. To do so, we binned the 
data into 60 s, 200 s, 600 s and 1000 s time bins. A visual inspection of the lightcurves did not 
reveal variability for any of the discrete sources at a level significantly exceeding the errors 
of the observation. 

Finally, we compared the X-ray point source population of NGC 3077 with optical B and 
a R c band images (see Sect. 2.2.1). Optical counterparts were not detected. An additional 
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SIMBAD query did not reveal known sources coinciding with any of the discrete X-ray 
sources. 



2.3. Diffuse Emission 

To obtain broad band images of diffuse emission, the discrete sources discussed in 
Sect. 2.2.2 were subtracted from the data over an area set by three times the correspond- 
ing PSF. The blanked areas were subsequently refilled by a Poisson distribution interpolated 
from nearby regions within the area of diffuse emission (task dmfilth in CIAO). The resulting 
images were then used to construct azimuthally averaged surface brightness profiles of each 
band, taking the infrared H-band peak as the center (a J20 oo = 10 h 03 m 19 s , 5 J20 oo = 68°44'02" 
derived from HST NICMOS data, Boker et al. 1999, see also Fig. 1 and 2b). Logarithmic 
representations of these profiles are displayed in Fig. 8. All bands are consistent with an 
exponential decline which is supported by corresponding fits to the data. Most of the X-ray 
photons are emitted in the medium band. The scalelength h of the surface brightness profile 
of this band (200±9pc) is thus similar to the one of the total band (h = 186±10pc). In con- 
trast, the soft energy band decreases more slowly (h = 296±49pc) and the hard band more 
rapidly (h = 135 ± 21 pc). The scalelengths quantitatively corroborate the morphological 
picture discussed in Sect. 2.1. 

With the help of the azimuthally averaged surface brightness profile of the total band, 
we were able to construct a radial volume density profile. For simplicity we assumed a spher- 
ical geometry of the hot coronal gas and decomposed the morphology into three-dimensional 
onion-like shells coincident with the two-dimensional annuli defined for the surface bright- 
ness profiles (for an example see the inset of Fig 9). The emission measure EM detected 
in each annulus k is the sum of the squared densities n within the corresponding shells i 
multiplied by their lines of sight l\ the general expression being 

M 

EM{k) = Y,(Uvr 1 n{i) 2 l{k) t . (3) 

i=k 

M denotes the outermost annulus, where the density can be easily calculated via n(M) = 
y/EM(M)/l(M) M . We further assume that the diffuse emission is emitted by a volume 
filling (f v = 1) hot collisional thermal RS plasma, with a temperature of 2 x 10 6 K (solar 
metallicity £ = 1 and solar element mixture), absorbed by a column density of 6 x 10 21 cm~ 2 
(for a justification see below). The relationship between the volume density of the hot 
gas and its measured emission measure and flux (see the surface brightness measurements 
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Sect. 2.3) was obtained by a simulation of the adopted plasma adjusted to the observational 
parameters and calibration. The volume density of the outermost shell is used as the starting 
value for an iteration towards the center of NGC3077; the l(k)i were re-calculated for each 
annulus k and shell %. Fig. 9 shows the resulting three dimensional radial volume density 
profile of the hot gas in NGC3077. The inset of this figure is an example for the method 
described using six shells. 

We fitted two models to the density profile (see also Sect. 2.4): an exponential (n(r) = 
ae~ r l h ) and a power law model (n(r) = (r/r c ) _,a ). The best fitting exponential has a 
scalelength of h = 271 ± 22 pc (Xred = 5 x 10~ 4 ), whereas the best power law has an index 
of = 0.57 ± 0.06 (x 2 red = 17 x 10~ 4 ). 

2.4. Evidence for Expanding Ha Shells Filled with Hot Gas 

A comparison of the raw photon maps with an Ha image shows that expanding Ha 
shells surround certain areas where diffuse X-ray photons are detected. This is also stressed 
by the adaptively smoothed representation shown in Fig. 2b. The situation is similar to 
what is expected from warm photoionized shells confining hot bubbles (Weaver et al. 1977). 
The morphological relationship of the images in the two different wavelengths is not due to 
spurious artifacts of the adaptive smoothing process. A convolution of the raw X-ray image 
with a fixed Gaussian kernel of ~ 10" is very similar to the results shown in Figs. 1 and 2. 

For a more quantitative analysis we constructed an intensity plot of azimuthally ar- 
ranged segments based on the unsmoothed X-ray data and an Ha map. In order to avoid 
confusion with bright H II regions, the inner 20" of NGC3077 were excluded from the anal- 
ysis. Some of the segments cover regions where the Ha morphology reveals expanding shells 
and some are in between. Fig. 10 supports the existence of Ha bounded hot bubbles: e.g., 
segments 7, 8, and 9 are more X-ray luminous than the surrounding segments 6 and 10/11 
where the Ha emission peaks. In contrast, the latter two segments cover regions where the 
rim of one particular expanding superb ubble (R2, see Fig. 11) is present. The same behavior 
is found for segment 5 where the X-ray intensity increases while the Ha flux goes down. 
Segments 12 to 3 to the south east of NGC3077 are strongly affected by absorption (see 
Sect. 2.1) which makes a similar comparison more difficult. The confinement of the hot gas 
is corroborated by the morphology of an H I feature to the west of NGC 3077 which resembles 
the rim of a supergiant shell and which appears to be a barrier for the hot, coronal gas (see 
Fig. 2c). 

We can also test this hypothesis by examining the density profile of the hot gas: For 
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a spherically-symmetric, outflowing wind, theoretical models predict n oc r _/3 with (5 — 2 
and T oc r -4 / 3 (n: volume density, r: distance along outflow, T: temperature; Chevalier 
& Clegg 1985). The same dependencies would be observed for a conical outflow with a 
constant opening angle (S97). A cylindrical geometry, however, results in a constant density 
all along the outflow ((3 = 0). As we show in Fig. 9, a power law with an index of (3 — 2 does 
not fit the volume density profile of the hot gas detected in NGC3077. As dwarf galaxies 
have low gravitational potentials, one naturally would expect that these systems are to first 
order spherically-symmetric, but even the best fitting power law ((3 = 0.57, Sect. 2.3) is 
a poorer fit compared to an exponential decline. In addition, we fitted MeKaL and RS 
thermal plasma models to a circular region 10" in radius centered on NGC 3077 and to three 
surrounding consecutive annuli with the same widths. The resulting temperatures varied by 
~ 30% around 2.1 x 10 6 K but no clear radial temperature gradient within the hot outflow 
is measured. This implies that it is unlikely that NGC 3077 currently forms a powerful hot 
galactic outflowing wind similar to what is witnessed, e.g., in M82 (well fit by a power law 
with an index of (3 = 0.9 and with varying opening angle Heckman, Armus, & Miley 1990; 
Bregman, Schulman, & Tomisaka 1995) or NGC 253 ((3 = 1.3). We conclude that the hot 
gas is indeed confined to giant Ha shells. 

Based on this analysis we consider it justified to divide the diffuse X-ray emission of 
NGC 3077 in regions (R2-R7) coinciding with expanding supergiant shells (see Fig. 11). One 
region (Rl) is defined to contain all of the X-ray emission within NGC 3077. We extracted 
ACIS-S3 spectra for all regions. Given their positions on the CCD, all spectra were calibrated 
using the appropriate RMFs and ARFs. Additionally, we defined background regions with 
the same sizes and on similar CCD rows. To check how significantly the calibration varies over 
the size of the aperture, we also constructed intensity-weighted response files. The differences 
to the spectra where the calibration files were not weighted were found to be negligible. The 
total X-ray emission of NGC 3077 (Rl) has enough counts to allow \ 2 statistics to work 
properly. We were therefore able to simply subtract the corresponding background counts 
from the source data. In contrast to Rl, the number of counts per spectral bin from regions 
R2 to R7 is much lower and we applied the Cash maximum likelihood algorithm accounting 
for Poisson statistics. As a consequence, we were not able to subtract the background from 
the data and we performed simultaneous fits of the background (power law = powlawld) 
and of the source plus background emissions ([photoelectric absorption x RS model] + 
power law). In addition, we fitted MeKaL thermal plasma models to the spectra (Mewe, 
Gronenschild, & van den Oord 1985; Kaastra 1992 and updates 4 with Fe-L calculations from 
Liedahl, Osterheld, & Goldstein 1995): ([photoelectric absorption x MeKaL model] + power 
law). The errors of the fits were again obtained by simulations as is described for the point 
sources in Sect. 2.2.2. The differences of the RS and the MeKaL models are considered to 
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provide a measure for systematic errors. The optimization method was the same as for the 
point sources (Monte-Powell) and again we adopted solar metallicity and element mixture. 
The resulting best fit parameters are tabulated in Table 3 (Nh'- absorbing column density, T: 
temperature, F x bs : absorbed flux, Fx' unabsorbed flux, Lx'- unabsorbed X-ray luminosity; 
errors are given at a 68% confidence level). The derived values for Nh are rather high while 
temperatures are low. Unfortunately, this has the effect that the conversion from absorbed 
to unabsorbed fluxes is achieved by rather large correction factors (~ 30 — 100). The derived 
unabsorbed X-ray fluxes and luminosities are therefore to be taken with caution. Note 
that plotting the hardness ratios of the different regions in a hardness ratio diagram (the 
shaded region in Fig. 14) reveals that the diffuse emission is very unlikely to be due to a faint 
population of underlying point sources with power law spectra. Thermal plasma models, 
however, are in good agreement with the observed hardness ratios. 

It is difficult to decide, which plasma model is a better fit to the spectra of the individual 
shells. Even for the best fits to Rl (the entire galaxy) goodness-of-fit indicators yield very 
similar values for both the RS and the MeKaL model. In Fig. 12c and d we show this 
spectrum with the best fits overlaid. Additionally, we calculated confidence regions for both 
models (Fig. 12b): the models agree well within a temperature range of T = 0.15 — 0.25 keV 
(equal to 1.7 — 2.9 x 10 6 K) and an absorbing column density range of Nh = 5 — 8 x 10 21 cm -2 
(boundaries of the 68.3% confidence levels). Similar plots are shown in the lower panels of 
Fig. 13 for the shells R2 to R7. 

The best fitting absorbing column densities are quite high for the entire galaxy compared 
to its measured H I column density (mean H I column density: ~ 2 x 10 21 cm~ 2 , peak: 
~ 4 x 10 21 cm~ 2 , W02). However, molecular gas is present in NGC3077 and can be detected 
by CO transitions at the center. In addition, the confidence regions of the best fit to Rl 
(Fig. 12b) shows that at least the MeKaL model constrains the lower limit for absorption very 
poorly and absorbing columns as low as say 10 21 cm~ 2 are well within the la uncertainty. 
The RS model, however, sets a lower 1 a limit of ~ 5 x 10 21 cm -2 which is difficult to reconcile 
with H I and CO data. 

From the normalization of the spectra (see footnote a in Table 2), we were able to derive 
the densities and the emission measures of each shell, assuming that the source geometry is 
spherical. It is also assumed, that the electron and the proton densities are equal, averaged 
over the entire volume (< n e >~< n p >). Note that the measured densities n rms must be 
converted to n e by using the appropriate filling factor /„ which is defined as r? rms = f v n 2 e . The 
pressures P/k of the shells can then be derived via P/k = 2n e T (k: Boltzman's constant, 
T: temperature). In addition, the hot gas mass Mhot = n e m p V (m p : proton mass, V: shell 
volume), the thermal energy E t h = 3n e VkT, the cooling timescale t coo i = E t h/Lx, the mass 
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deposition rate M coo i = Mhot/t CO oi, and the mean particle velocity < v\ lo t >= \/2E th / M hot 
are listed in Table 4 which also shows their dependence on /„ and on the metallicity of 
the hot gas (£, in solar metallicity units). The equivalent diameters d eq were derived as 
being the diameters of circles which have the same areas as the regions under consideration 
(d eq = a/4 Area/n). The mean lines of sight are obtained by 2/3xd eq (cylindrical projection). 
Errors imposed on the physical gas parameters by the source geometry are estimated to be 
of order 30%. 

For the individual shells we derive temperatures of ~ 1.3 to 4.9 x 10 6 K, hot gas masses 
of order ~ 10 5 £~ a5 /°- 5 M , and particle volume densities of ~ 0.05 — 0.9 £~ a5 /~ a5 cm -3 . 
The pressures of the hot gas derived for the individual shells in NGC3077 are P/k m 2 — 
30 x 10 5 £ _a5 /~ a5 Kcm _3 . This is in contrast to typical pressures of the ISM in the Milky 
Way which are of order a few thousand (e.g., Wolfire et al. 1995). It is this overpressure of 
the hot gas which drives the expansion of the superbubbles. 

Taking the hot gas masses M hot and the age of the shells t she u (Sect. 3.2.1 and Ta- 
ble 4), we can estimate the rate of hot gas deposed to the halo. We derive M hot /t she u ps 
0.1 - 0.5^- - 5 f°- 5 MQyT- 1 which is similar to the warm, ionized gas (M warm /t s h e ii = 0.6 x 
(fv,warm/0.l) - 5 M yr -1 , M98) and a few times greater than the current star formation rate 
(~ 0.06 Meyr" 1 , W02; value adapted to a distance of NGC3077 of 3.6 Mpc). 

The value of the filling factor /„ is uncertain. Hydrodynamic simulations and Chandra 
observations of superwind galaxies (Strickland & Stevens 2000; Strickland et al. 2002) show 
that /„ is about 0.1-0.3. The hot gas in NGC3077, however, is still confined to giant shells 
and the X-ray emission emerging from these regions does not appear to be brighter in the 
vicinity of the Ha emitting rims as compared to the interior of the shells (see above). A more 
quantitative criterion is that the time t c needed to cross a superbubble with the sound speed 
c of the plasma is less than the age of the shell. The sound speed is given by c = ^j2^kT /m p 
where 7 = 5/3 (S97). Taking a mean temperature of 2 x 10 6 K and a diameter of the bubble 
of 300 pc (see Tables 3 and 4) we derive c 240 km s" 1 and t c 1 Myr. t c is a factor 2-10 
less than the ages of the bubbles and the hot plasma had enough time to thermalize within 
the cavities. A volume factor close to unity therefore appears to be more satisfying than the 
lower ones detected in outflows. 

There are some differences in the properties of the individual shells. R7 and R4 were 
fitted to have high absorbing column densities and low temperatures relative to the others. 
The combination of these parameters lead to relatively high correction factors for absorption 
and therefore high unabsorbed fluxes and luminosities. The opposite is the case for R3 
and R6, the best fits suggest lower absorbing column densities and higher temperatures. 
The highest number of photon counts is observed in R2 which corresponds to the highest 
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absorbed flux for all shells. However, after correcting for absorption, the unabsorbed X-ray 
luminosity is in the mid range of the shell luminosities. In addition, R7 and R4 have the 
highest normalizations which results in relatively high volume densities, masses, thermal 
energies and pressures of the hot gas. Again, R6, R3 and also R5 show relatively low values 
for these parameters. The reason for this might be that R7 and R4 are located very close 
to the stellar disk of NGC 3077 (see Fig. 2b) which is centered on R3 with a position angle 
of ~ 45° (taken from the LEDA catalog Prugniel & Heraudeau 1998). In this case, R3 
itself, which seems to be absorbed by a layer of rather low column density, must be located 
closer to the observer as compared to R4 and R7, in particular in front of the extended 
H I feature in the south-east and the molecular gas in the central region of NGC 3077 (see 
Fig. 2c). Alternatively, NGC 3077 maybe inclined in such a way that the far side of the disk of 
NGC 3077 points towards the south-east. However, the distribution of H I within NGC 3077 
is rather chaotic and is not in the form of a regular disk along with its optical counterpart. 
Significant contamination with non-thermal X-ray photons is not supported by the spectral 
data of R3 (and all the other shells); above an energy of ~ 1.5 keV no emission exceeding 
the background can be detected (see also the shaded region in the hardness ratio diagram 
Fig. 14). 



2.5. The Total Spectrum 

We extracted a total X-ray spectrum of NGC 3077 including both point sources and 
diffuse emission (Fig. 15). In order to verify the results obtained for the individual spectral 
fits, we added the following best fit model components and compared the results to our 
measurement. Point sources: RS plasmas for SI, S5, and S6, power law for S2 and S3, 
black body for S4 (see Sect. 3.1 and Table 2); Diffuse Emission: RS/MeKaL plasma for Rl 
(see Tab. 3). For all models, we adopted the absorbing column densities as derived by the 
corresponding fit. For energies < 0.8 keV and > 1.2 keV the total model provides a good fit 
to the data. In the 0.8keV< E < 1.2 keV energy range, the total model shows a count rate 
which is ~ 5 x 10~ 3 cts s _1 keV -1 too high. This slight offset is most likely due to the fits to 
the point sources, where we modeled the background by a power law. In Fig. 15 we show the 
added best fitting RS/MeKaL plasma models for the regions R2 to R7 as well. The derived 
combined flux of the regions R2 to R7 is a factor of ~ 1.6 lower than the diffuse X-ray flux 
of the entire galaxy (Rl) in the case of both, the RS and the MeKaL models. The total 
unabsorbed luminosity of NGC 3077 is 2 — 5 x 10 39 erg s _1 , depending on the plasma model 
used. 
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2.6. Is there an X-ray Population in the Tidal Tail close to NGC 3077? 

NGC3077 finds itself in interaction with M81 and M82 (the M81 triplet of galaxies). 
Large H I tidal tails cover the region around these galaxies (see, e.g., Yun, Ho, & Lo 1994). 
To the east of NGC 3077, a massive tidal H I complex is present and molecular gas and star 
formation as traced by CO and H II regions have been found in this tidal feature (e.g., see 
Walter & Heithausen 1999, W02, Karachentsev, Karachentseva, & Boerngen 1985). The roll 
angle and the aimpoint of the Chandra observations presented here were carefully chosen 
in a way such that in addition to NGC 3077 the bulk of the H I emission and most of the 
extragalactic H II regions are within the ACIS-S3 field of view (see Fig. 5). 

In addition to the point sources associated with NGC 3077, we detect 38 point sources in 
the ACIS-S3 field of view, 24 of which were detected by both the wavdetect and the celldetect 
algorithm (see Sect. 2.2.1). We compared the source list (Table 1) to the SIMBAD database 
and found that three sources correspond to known stars: two form the bright double star 
HD 8667 to the north and one source is coincident with a star detected by the Hubble Space 
Telescope pointing at the halo of NGC 3077 (Sakai & Madore 2001). The counts for each 
X-ray point source and band are listed in Table 1. 

Fig. 14 shows a hardness ratio plot of all point sources in the field of view which were 
detected by both source detection algorithms. Most of the sources are consistent with mod- 
eled power law spectra. Clearly offset are the three detected stars as well as S4, S5, and 
S6 within NGC 3077, which are substantially softer (see Sect. 3.1.1 and 3.1.3). This implies 
that we find no evidence for a bright population of soft sources in the tidal tail (e.g., SNRs). 

To check whether the detected X-ray sources belong to the M81 triplet or if they 
are background objects we compared the data to the results from the Chandra Deep Field 
South project (CDFS; Giacconi et al. 2001, in the following G01). To do so, we re-reduced 
our Chandra data in the same manner as G01. Within two new energy bands (G01 soft: 
0.5keV< E <2.0keV, G01 hard: 2.0keV< E <7.0keV), the data were binned to a pixel size 
of 1" and were corrected for vignetting effects. Fluxes were calculated regarding the Galactic 
foreground column density along the line of sight to NGC 3077. For a full description of the 
data reduction process see G01. 

Excluding the identified stars and the sources within NGC 3077, we constructed log N(> 
S) — log S plots (N: cumulative number of sources brighter than S per square degree, S: 
flux in erg s _1 cm~ 2 ) for both of the G01 bands. In Fig. 16 we compare our results to those 
obtained by G01 for the CDFS. All data are well described by the G01 relations. We do 
not observe a significant overdensity of X-ray sources within the tidal tail around NGC 3077 
down to our detection limits. This is in agreement with scaling the number of the point 
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sources within NGC 3077 by the ratio of blue and Ha luminosities of NGC 3077 and the 
"Garland" region (see Karachentsev, Karachentseva, & Boerngen 1985) which occupies a 
considerable fraction of the tidal tail {L% arland / L% GCm77 = 5.2 x 1O 6 L B0 /1.1 x 1O 9 L B0 « 
1/200; Fg% rland /FgGC3077 = Lg8 x iQ-^erg s" 1 cm- 2 /3.75 x 10- 12 erg s" 1 cm- 2 « 1/20; 
Sharina 1991, de Vaucouleurs et al. 1991, Walter, Martin, & Ott 2003). Assuming that 
these ratios are a measure for the ratio of the number of X-ray point sources, we expect 
to detect no X-ray point source in the Garland region. This estimate is in agreement with 
our analysis and in contrast to the results of Bi, Arp, & Zimmermann (1994). They find 
an overpopulation of X-ray sources in the tidal tail of a factor of two, based on their 7ks 
ROSAT data of a field centered on NGC 3077 and 10' in size (roughly corresponding to the 
ACIS-S3 field of view). 



3. Discussion 

3.1. The Nature of the Discrete X-ray Sources in NGC 3077 

3.1.1. Supernova Remnants 

The spectrum of SI (Fig. 6) shows identifiable line complexes. SI is spatially coincident 
with radio continuum emission, a CO peak and a strong absorption feature in the optical. 
The Paa (Boker et al. 1999) and the radio continuum morphology are similar, which suggests 
that the radiation arises from the same region and is due to shocked gas. This would agree 
with the fact that we derive a radio continuum spectral index of a — 0.48 (def.: S oc z/- Q , 
data from Niklas et al. 1995; Condon 1987) which is consistent with spectral indices found 
by Hendrick & Reynolds (2001) for a sample of SNRs in the Large Magellanic Cloud. 

At the distance of NGC 3077, the slightly extended source S5 has a size of 19pcx23pc. 
The hardness ratio plot (Fig. 14) reveals, that S5 and S6 are rather soft sources compared to 
background objects and to power law models. From the fitting process we cannot reliably 
decide which of the three models (RS, PL, BB) provides the best fit to the data of SI, S5, 
and S6. However, the photon indices of the PL model are too high to be explained by any 
known celestial source. In addition, the alternative BB models yield temperatures which are 
in excess of at least a factor of two to models of white dwarf or isolated neutron stars (see 
Kahabka & van den Heuvel 1997; Kahabka 1999; Haberl, Pietsch, & Motch 1999; Burwitz et 
al. 2001; Zampieri et al. 2001). For these reasons, we conclude that SI, S5, and S6 may be 
(young) SNRs and emit X-rays from a thermal plasma. They are close to the actively star 
forming regions in NGC 3077 and to the molecular complexes (see Fig. 2). From the best 
fitting RS and MeKaL plasma models we derive the parameters given in Tables 2 and 5. 
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3.1.2. Accreting Objects 

S2 and S3 exhibit a somewhat similar spectrum with no clear energy peak. The best 
fitting RS and BB spectra result in high temperatures. PL fits are supported by the hardness 
ratio plot (Fig. 14) where the hardness ratios of S2 and S3 are similar to those derived for 
background objects and are consistent with simulated PL models. The fitted absorbing 
column densities differ by a factor of ~ 3; it is likely that they are unrelated objects, in 
spite of their small angular separation (1"). Due to the low S/N, we can only speculate that 
these sources may be accreting objects. Judging from their luminosities, the nature of these 
objects can be both, low or high mass X-ray binaries. Both sources might be background 
AGNs, too. Especially S2 might be such an object, due to the steepness of its spectral 
index (7 = 1.65; cf. sources in the CDFS: 7 = 1.70 ± 0.12, see G01) and its high absorbing 
column density which might place S2 behind NGC3077. Statistically, given the density of 
background sources in a typical field (e.g., the CDFS), a background object coinciding with 
the diffuse emission of NGC 3077 is possible. 

3.1.3. The Supersoft Source 

The point source with the softest spectrum is S4. It is located on the rim of an Ha 
feature. Its lack of line features and emission harder than 0.8 keV suggests that the spectrum 
is best described by a black body model (T = 9.4 x 10 5 K equivalent to 81 eV). A power law 
index of 5 would be too steep for any Galactic or extragalactic source and the best RS plasma 
fit features the Fe-L line complex near 0.9 keV which is not observed. The temperature and 
the X-ray luminosity of 9.2 x 10 36 erg s _1 of the BB spectrum puts this object within the 
range found for other supersoft sources (Kahabka & van den Heuvel 1997; Swartz et al. 
2002). 

It is possible to derive the radius of an object from the amplitude of the BB spectrum 
(see footnote b in Table 2). For the distance of NGC 3077 we find a radius of R ~ 1600km. 
This is much lower than the radius of a typical supersoft source, a hydrogen burning white 
dwarf (WD; R ~ 6000 — 9000 km; Kahabka & van den Heuvel 1997). However, since the 
hydrogen burning phase is confined to a small region on the surface of the WD, the derived 
value might be in accordance with such an object. 

An alternative interpretation would be that S4 is an isolated neutron star. The radius 
of a typical neutron star is about 20 km. If we adopt this typical radius, then S4 would 
have a distance of only 45kpc, which puts it in the halo of the Milky Way. Located at this 
distance, its luminosity would drop to Lx = 1-4 x 10 33 erg s -1 which is rather high for an 



isolated neutron star but not impossible (see, e.g., Zampieri et al. 2001). The best fit for the 
column density is similar to the Galactic foreground hydrogen content and is even consistent 
with a value of zero within the errors (Fig. 7). 

One can estimate the expected number of supersoft sources by a comparison of NGC 3077 
with M81, its neighbor galaxy. M81 was observed with Chandra's ACIS-S3 CCD for about 
50 ks. This observational setup and integration time is comparable to the X-ray observa- 
tions discussed here. Based on that data, Swartz et al. (2002) find 9 supersoft sources in 
M81. The absolute blue magnitudes (M81: -19.8 ±0.14 mag, NGC3077: -17.1 ±0.14mag; 
de Vaucouleurs et al. 1991) show that M81 is about a factor of 10 more luminous than 
NGC 3077. Scaling the number of supersoft sources by this factor, one would expect one 
supersoft source in NGC 3077 which is what is observed. 

3.2. The Diffuse X-ray Emission of NGC 3077 

3.2.1. Pressure-Driven Bubbles 

Theoretical models of expanding supershells are generalized models of winds from mas- 
sive stars streaming into the ambient ISM (see, e.g., Castor, Weaver, & McCray 1975; Weaver 
et al. 1977; Mac Low & McCray 1988). In the following, we will use the set of equations 
derived by Chu & Mac Low (1990) who used the EINSTEIN band (0.2 - 4.0 keV) over which 
to calculate the X-ray emissivity of hot gas. For temperatures of about 2 x 10 6 K, no substan- 
tial emission above 4.0 keV is expected, so their results should be applicable to the Chandra 
data discussed here. The theoretical X-ray luminosity of a pressure-driven superbubble is 
given by 

L X = 3.29 x 10-erg { /(r) t%» (4) 

where £ is the metallicity of the hot gas with respect to solar, L mec h^ the mechanical 
luminosity in units of 10 37 erg s -1 , n am b the number density of the ambient material in cm~ 3 , 
and t she iifi the lifetime of the starburst in Myr. The integral I(t) = (125/33) — 5-r 1 / 2 + 
(5/3) r 3 — (5/11) r 11 / 3 only depends on the dimensionless temperature r = T min /T center given 
a cutoff temperature T min and the temperature at the center of the shell T center . Expressed 
in observables: 



r _ n i fi r-8/35 -2/35,6/35 
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The ambient density n am b is related to the electron density in the interior of the shell n e by 

n e kT 

riamb ~ 5— (6) 

where k is Boltzman's constant, T the temperature, m p the proton mass, and v exp the 
expansion velocity of the shell. 

The density-dependent mechanical luminosities L mec h/n am b, the ages t s h e ii, and the 
velocities v v J° v ^ of the shells were previously derived by M98 using long-slit echelle spectra 
(see Table 4). Adopting solar metallicity and the best fit values from Table 3, J(r) is about 
2.8 and the resulting predicted densities of the ambient material, the mechanical luminosities 
Lmech, and the predicted X-ray luminosities are listed in Table 4. Note, that shells A and J 
in M98 have a somewhat larger size compared to R3 and R7. However, this should not have 
any significant impact on the predicted numbers. 

A comparison of the theoretical values with the ones directly derived from the Chandra 
data (Tables 3 and 4) reveals that the observed luminosities are in reasonable agreement with 
the basic theory developed by Weaver et al. (1977). In their model, "boiled-off" material 
from the rim of the superbubble is the only mechanism for mass-loading. Therefore, the hot 
gas which fills the cavities is expected to originate only from the evaporation process of the 
rim as well as from the stellar ejecta themselves. We estimate the masses of both processes. 
Weaver et al. (1977) derive the evaporated mass rate dMb/dt to 

d _Ml = c <T >5 / 2 _^!_ (7) 
dt R-r w 

where C = 4.13 x 10~ 14 (cgs units) is a material-dependent constant, < T > is the 
mean temperature of the hot gas, R the radius of the superbubble, and r the radius of a 
sphere where the stellar ejecta has not thermalized yet. Given the large sizes, ages, and the 
plasma sound speed of the superbubbles in NGC3077 (see Sect. 2.4 and Table 4) we expect 
that all the material is already thermalized, i.e., r = 0. By using the derived values given 
in Table 4 and integrating over the age of the individual bubbles we estimate the amount of 
the evaporated mass to ~ 20 — 60% for R3, R5, and R6 and ~ 1 — 10% for R2 and R7 of 
the hot gas masses. The uncertainties are given by the different plasma models. 

In addition, we used the STARBURST 99 models provided by Leitherer et al. (1999) 
to estimate the injected mass from the combined effect of stellar winds and SNe. To do so, 
we estimate the star formation rate (SFR) within each superbubble by scaling the measured 
mechanical luminosity (Tab. 4) by the asymptotic values of the STARBURST 99 models with 
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continuous star formation (Salpeter IMF, upper mass cutoff: f00M Q , lower mass cutoff: 
1M , solar metallicity). The SFRs over the lifetimes of the superb ubbles are derived to 
0.02 — 0.06 M Q yr^" 1 for each superbubble except for R7 (~ 0.4 M yr -1 ). These values are 
subsequently used to scale the mass loss predicted by the theoretical models integrated over 
the ages of the individual superbubbles. As a result we estimate that about 30 — 50% of the 
hot gas within R3, R5, R6, and R7, and about f — 10% within R2 may originate directly 
from stellar ejecta. 

A mass balance shows that additional mass-loading besides evaporation of the rims and 
stellar ejecta, e.g., by evaporation of halo or disk clouds, is definitely needed for the largest 
superbubble R2 in order to account for > 80% of its hot gas mass. R7 needs some ~ 40% 
of additional mass-loading and the other superbubbles R3, R5, and R6 can be explained 
entirely without this mechanism. 

Adding up the estimated SFRs of the superbubbles obtained by the STARBURST99 
models yields ~ 0.6 M yr _1 which is about an order of magnitude higher than the current 
SFR (0.06 M Q yr" 1 , W02, based on Ha and Paa flux measurements). The hot and warm gas 
deposition rates lead to similar discrepancies with the current SFR (Sect. 2.4). Note that 
the 1.4 GHz radio continuum flux of W02's VLA data (~ 42 mJy) results in a similar current 
SFR of ~ 0.08M Q yr _1 (conversion factor taken from Haarsma, Partridge, Windhorst, & 
Richards 2000). Also, the far-infrared luminosity of NGC3077 (3.2 x 10 s L , Yun 1999, 
adapted to a distance of 3.6 Mpc) can be converted to a SFR of ~ 0.06 MQyr^ 1 (conversion 
factor taken from Kennicutt 1998). The differences of the STARBURST99 and current SFR 
estimates might indicate that the star forming activity of NGC 3077 has decreased with time 
over the last few Myr (for the ages of the shells, see Table 4). 



3.2.2. Impact of Superbubbles on the Intragroup Medium 

The superbubble shells will accelerate at radii where the ambient gas density drops 
rapidly. As hydrodynamic instabilities break up the shells, the hot gas will rush outwards 
in a galactic wind. An obvious question is whether the outflowing material will return to 
NGC 3077 or not. The answer depends on several factors. A necessary condition is that 
the wind has enough energy to escape from the gravitational potential. The wind can still 
stall, however, if the hot wind cools and/or does a great deal of work pushing an extended 
gaseous, galactic halo aside. The cooling times of the shells (2-20 Myr, see Table 4) exceed 
the flow time (1-10 Myr) so radiative losses do not direct the dynamical evolution at this 
stage. Hence, we discuss only the impact of the gravitational potential and gaseous halo. 
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The tidal interaction with M81 and M82 (e.g., Yun 1999) raises questions about the 
reality of mass models for NGC3077. M98 assume a H I rotation speed of 40 km s _1 across 
the system. However, higher angular resolution H I data (W02) subsequently showed the H I 
velocity spread within the optical body of NGC 3077 to be 90 km s _1 (la velocity dispersion 
of ~ 25kms- 1 ). W02 detect 6.1 x 1O 8 M of H I in the tidal tail (including NGC 3077) 
but only 1.3 x 1O 8 M of H I in NGC 3077 (values adjusted to a distance of 3.6 Mpc), and 
no rotation was detected in the latter. The stellar mass of NGC 3077 is estimated to be 
1.4 x 1O 9 M based on the blue luminosity and M/L B = 1 (de Vaucouleurs et al. 1991). 
Hence, the total baryonic mass is probably 1.5 — 2 x 10 9 M . 

Within a simple model, we scale the dark mass M DM of NGC 3077 by its visible mass 
M vis according to the empirically derived relation (Persic, Salucci, & Stel 1996; Silich & 
Tenorio-Tagle 2001) 

M DM = 3.74 x 10 8 (y^) 071 M (8) 

and estimate M^m ~ 1O 1O M leading to M^M/M vis ~ 10 (in agreement with the 
estimate of Brouillet et al. 1991, based on M/L ps 9). The complicated history of NGC 3077, 
which might have had a much larger (or smaller) mass in the past, however, amplifies the 
uncertainty in the ratio of dark to baryonic matter. 

Given an isothermal dark matter halo, the escape velocity is a function of radius r and 
rotation speed v rot . It is described by 



Vesc = V rot (l + In j 



(9) 



with a maximal halo radius r max of 



W = °- 016 (^) ^ 2/3k P c ( 10 ) 

(h: dimensionless Hubble constant, adopted to h = 0.7 Mac Low & Ferrara 1999; 
Ferrara & Tolstoy 2000; Silich & Tenorio-Tagle 2001). For NGC 3077 we derive an escape 
velocity of ~ 110 km s -1 at the center and ~ 80 km s _1 at lkpc above the disk. Expressed 
as an escape temperature, the bubble gas needs to exceed T = 2 — 9 x 10 5 K to be assured of 
escape. This condition is clearly met by the hot bubbles (T « 1 — 3 x 10 6 K). In contrary, the 
expansion speed of the warm shells is about equal to the maximum escape velocity predicted. 
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As a second condition, the expanding bubbles must provide the work to push the halo 
gas aside. The velocity of a shell v s is a function of the pressure of the hot gas Phot and 
the density of the ambient medium p am b'- v s = (Phot/ Pamb) ' 5 (Koo & McKee 1990; Silich & 
Tenorio-Tagle 2001). For decreasing density, e.g., a stratified atmosphere, the velocity of a 
shell increases as long as the driving pressure is constant. For an ultimate loss of material 
the density has to decrease sharply before the hot gas in the interior of the shell cools down. 
Koo & McKee (1992) show that for an exponential density profile the shell has to reach 
three times the exponential scalelength before it re-accelerates, fragments, and ultimately 
blows the currently stored hot gas out of the galaxy's gravitational potential. Assuming a 
Gaussian stratification, this effect will start at about one Gaussian scalelength. 

The geometry of the H I data is more complicated, but we argue that escape is clearly 
easiest where the gradient in halo gas density is steepest, which is roughly along the north- 
south axis (see Fig. 5). This is also the direction where the largest elongation of the diffuse 
X-ray emission is detected (Fig. 1). Towards the south, NGC3077 is connected to a tidal 
tail which implicates that the H I column density is very extended (~ 5kpc) at a constant 
Nhi ~ 1-5 x 10 21 cm~ 2 (Fig. 17). Towards the north, however, a continuous decline is 
observed. Hot gas streaming to the south will most likely not be able to overcome the large 
ISM pressure "barrier"; the cooling times of the hot gas within the shells are much lower 
than the time needed to reach the outer boundary of the southern tidal tail (see Table. 4). 
This is in agreement with the morphology of the diffuse X-ray emission (Fig. 1) which is 
considerably less extended towards the south compared to the north. 

We fitted a Gaussian and an exponential profile to the northern half of the H I col- 
umn density profile of NGC3077 and derive the respective scalelengths to h gauss ~ 0.5 kpc 
and h exp ~ 0.6 kpc (for comparison: the scalelength of the total diffuse X-ray emission is 
~ 0.2 kpc, Sect. 2.3). Both profiles provide an acceptable fit. Via N HI = \f 7 hx h gauss n 
(Gaussian) and Nhi = h exp no (exponential), we can derive the midplane volume density no 
of NGC3077. Assuming that the mean column density is Nhi ~ 2 x 10 21 , we consistently 
obtain n = 0.5 — 1.0 cm -3 for both models. This is in good agreement with what we predict 
based on Eq. 6 for the ambient densities of the shells (see Table 4). If the Ha shells do not 
decelerate significantly (current maximum velocity: v ~ 100km s _1 ), they will reach the H I 
scalelengths in ~ 5Myr just before substantial cooling takes place (Table 4). In addition, 
current star formation provides additional energy to the shells. It is likely that the shells 
and the hot gas indeed reach the northern H I scalelengths while maintaining supersonic 
speed. Eventually, the shells will fragment due to Rayleigh-Taylor instabilities and release 
their hot gas to the intragroup medium of the M 81 group of galaxies. 

Recently, Karachentsev et al. (2002) estimated the mass of the M81 group to 3 — 16 x 
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10 11 M , its rotational velocity to ~ 250 km s" 1 and its size to ~ 500 kpc. Assuming M81 to 
be the center of the group (distance to NGC3077 ~ 50 kpc), we estimate the escape velocity 
of the group at the position of NGC3077 to ~ 500 — 600 km s _1 . This converts into an 
escape temperature of ~ 2 — 3 x 10 7 K, about an order of magnitude higher than the current 
temperature of the hot gas: the hot gas of NGC 3077 is not able to escape from the M 81 
group of galaxies. 

In summary, we conclude that some of the hot gas will likely enrich the intragroup 
medium of the M 81 group to the north but not to the south of NGC 3077. The IGM beyond 
the M81 group, however, will not be affected by the developing superwind of NGC 3077. 

3.3. Comparison to M81 and M 82 

NGC 3077 is interacting with its neighbors, the spiral galaxy M81 and the starburst 
galaxy M82. A compilation of some general and X-ray properties of these three galaxies 
can be found in Table 6 (the ROSAT HRI X-ray data of M 81 has been taken from Immler 
& Wang 2001, the ROSAT PSPC data of M82 from S97). In contrast to NGC 3077, both 
of its neighbors harbor a strong point-like nuclear source, each of which is believed to be 
associated with a black hole system (e.g., Iyomoto & Makishima 2001; Kaaret et al. 2001). 
To make a proper comparison of the X-ray properties, we exclude these sources in the 
following discussion (resulting in columns 5 and 8 of Table 6). Furthermore, Immler & Wang 
(2001) model the diffuse emission of M81 by two thermal plasmas and a luminous power 
law component. The latter is supposed to comprise faint point sources, rather than truly 
diffuse gas, which is why we add this component to the point source flux and extrapolate 
the luminosities accordingly (column 6). 

The dominating X-ray emission emerging from M81 is in the form of point sources 
(91%). This is not observed in M82 and NGC 3077 where most of the X-ray photons are 
radiated by hot outflowing gas (95% and 85% respectively). The diffuse X-ray luminosity 
of NGC 3077 is more luminous than the one in M81 even in absolute terms (Table 6). As 
expected, the absolute X-ray luminosity of the point sources is higher for more massive 
galaxies. 

S97 discuss the physical properties of the galactic wind of M82. They find a total hot 
gas mass of M hot = 1.3 x 10 8 /°- 5 M o , a thermal energy of E th = 3.6 x 10 56 /°- 5 erg, and 
a mass deposition rate of M cool = O^Moyr" 1 (NGC 3077: M hot « 3 x 10 6 £~°- 5 jf 5 M , 
E th « 3 x 10 54 £~°- 5 ff 5 erg, M cooi « 0.1 M yr" 1 , see Table 4). The cooling times t cool of the 
coronal gas in M 82 are derived to be about 600 f®' 5 Myr, and the temperatures are some 
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5 x 10 6 K (NGC3077: t cool « 40 £~ - 5 jf 5 Myr, T « 2 x 10 6 K). Point source subtraction in 
the S97 data, however, is a difficult task, given the low ROSAT resolution. Compared to 
M 82, the starburst of NGC 3077 is weaker in absolute terms. The temperature of the hot 
gas in NGC 3077, however, is only a factor of two lower as compared to that in M 82 and 
very similar to the dominating hot plasma component of M81 (~ 1.7 x 10 6 K). In addition, 
the pressure of the hot gas in M 82 is even a factor of ~ 2 lower as compared to NGC 3077. 
This is an effect of the ~ 5 times lower density within the outflowing superwind in M 82. 

Normalized to the H I contents or the blue luminosities of NGC 3077 and M82 (which 
are potential measures for the "fuel" of the starbursts and the total masses of their host 
galaxies), the strong galactic superwind in M82 still carries a higher mass and thermal 
energy of the hot gas. The situation is different, however, when normalizing these quantities 
to the current SFRs. The SFR obtained from Ha measurements, the far-infrared luminosity, 
the thermal energy of the hot gas and its mass, all these parameters are some two orders 
of magnitudes lower in NGC 3077 as compared to M 82. This result is somewhat intriguing 
given the very different physical states of the strong galactic superwind in M 82 and the hot 
gas confined in superbubbles in NGC 3077. Note that the star formation efficiency based on 
molecular gas (SFR/M H2 ) of NGC 3077 and M82 are approximately the same, too (W02). 



4. Summary and Conclusions 

In this study, we present deep (53 ks) Chandra observations of the dwarf starburst galaxy 
NGC 3077 and its environment. Our results are the following: 

1. The main X-ray emission emerges from a diffuse soft component lkpc in size. In 
addition, 6 point sources are detected. The point sources contribute 8-21% to the 
total X-ray luminosity, the diffuse emission 79-92% (depending on the plasma model). 

2. The diffuse X-ray emission at the center of NGC 3077 is hardened due to photoelectric 
absorption of cooler gas as traced by high-resolution H I and CO observations. The 
diffuse emission is well described by collisional hot thermal plasma models. 

3. A comparison of the diffuse X-ray emission with the Ha morphology of NGC 3077 
reveals that regions containing hot gas are confined by expanding Ha shells. The 
confinement is corroborated by the volume density profile which cannot be explained 
by a spherically-symmetric freely expanding wind. 

4. For the individual superbubbles we derive the temperatures of the hot gas to ~ 2 x 

10 6 K, the masses to ~ 10 5 £ -0 - 5 /°- 5 M , the particle volume densities to ~ 0.1 £ -0 - 5 /~ a5 cm -3 , 



-25 - 



and the luminosities to ~ 10 38 erg s" 1 , in agreement with theoretical models. The mass 
of the stellar ejecta and of material evaporated from the rim of the shells balances the 
mass of the hot gas within four shells. For two superbubbles, including that with the 
highest X-ray brightness, other mechanisms for mass-loading must be invoked. The 
high pressures (P/ k ~ 10 5-6 £"°' 5 f^ ' 5 K cm -3 ) of the hot gas are believed to drive the 
expansion of the Ha shells. The entire diffuse X-ray luminosity is ~ 3.0 x 10 39 erg s" 1 ; 
~ 1/3 of the flux cannot be attributed to individual shells. 

5. The starburst properties of M82 and NGC3077 are quite different. In contrast to the 
superwind observed in M 82, the hot gas in NGC 3077 is still confined by the supershells. 
Also the physical parameters of the hot gas such thermal energy and mass of the hot 
gas as well as the current star formation rate are much higher in M 82 in absolute terms 
(about two orders of magnitude). The temperature of the thermal plasma, however, 
is remarkably similar in M82, NGC 3077, and in the non-starburst galaxy M81. The 
mass and thermal energy of the hot gas in M 82 is still an order of magnitude higher 
compared to NGC 3077 when normalizing these quantities to the galaxies' H I masses 
or blue luminosities. Relative to their current star formation rates, however, both 
starbursts carry about the same amount of hot gas mass and thermal energy despite 
of the very different morphological and physical states of their starbursts. 

6. Spectral analysis of the six point sources reveal that NGC 3077 hosts three young 
supernova remnants, two X-ray binaries, and one supersoft source. The first five 
sources are located close to the center of NGC 3077, the latter one in the halo. 

7. The X-ray source population in the tidal arm close to NGC 3077 is consistent with the 
extragalactic background population of the Chandra Deep Field South. An additional 
population of X-ray binaries or SNRs with fluxes >6x 10~ 16 erg s" 1 cm" 2 (0.5-2 keV) 
can be excluded. 

8. Hot gas will likely escape from the gravitational potential of NGC 3077 towards the 
north once the superbubbles reach sizes larger than the scalelength of the neutral gas 
(~ 0.5 kpc). To the south, however, the hot gas may not be able to overcome the 
pressure of the environmental, extended ISM and will return to its host. The gas 
driven out to the north will be stored in the intragroup medium of the M 81 group. 
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Fig. 1. — X-ray emission of NGC3077 in different energy bands (logarithmic scale). The 
data shown here are adaptively smoothed. Soft: 0.3 keV^ E ^0.7keV, Medium: 0.7keV^ 
E <l.lkeV, Hard: l.lkeV< E ^6.0keV, and Total: 0.3 keV^ E ^6.0 keV. The locations 
and definitions of the point sources are labeled in the total band image and the white plus 
refers to the NICMOS H-band peak (see also Fig. 2). 
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Fig. 2. — a) Three-color composite of the X-ray bands soft (red), medium (green), and 
hard (blue), b) Three-color composite of Ha (red), the total point source subtracted X- 
ray emission (green), and the H-band image (blue) of NGC3077. All images have been 
logarithmically scaled. Note that Ha shells are filled with hot gas producing X-ray emission. 
This image forms the basis of the definition of the different regions in Fig. 11. c) Three-color 
composite of the atomic hydrogen (H I, red), the diffuse X-ray emission (green), and the 
CO content (blue). The X-ray emitting region is wedged between the cooler phases of the 
interstellar medium. 1' corresponds to 1 kpc. d) The smoothing scales computed by the 
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adaptive smoothing algorithm which were used to convolve the diffuse X-ray emission. All 
panels are on the same spatial scale. 
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Fig. 3. — Hardness ratio images of NGC3077 using the normalized adaptively smoothed X- 
ray images. The ratios are defined as HRl=(Soft— Medium— Hard)/(Soft+Medium+Hard) 
(top) and HR2=(Soft+Medium-Hard)/(Soft+Medium+Hard) (bottom). Both images 
show the same region. The contours change their color at HR=0 and they are separated by 
0.1. Point sources are labeled according to Fig. 1. 
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Fig. 4. — Grey scale: High-resolution (6'/5) H I column density map of NGC 3077 in units 
of 10 21 cm -2 . Black contours: The X-ray hardness ratio image HR1 (same contour levels 
as in Fig. 3). Grey contours: CO (1-0) line emission (CO data taken from W02). 
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Fig. 5. — The raw ACIS-S3 data superposed with H I contours of the tidal tail. The contours 
start at a column density of 2.76 x 10 20 cm~ 2 and are equally spaced by the same amount 
(resolution: 40"). Point sources detected by the wavelet algorithm (wavdetect) are marked 
with small ellipses. The numbers correspond to those in Table 1. Additionally, we show the 
definition of the diffuse emitting regions evaluated in Table 3. 
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Fig. 6. — The ACIS-S3 spectrum of the point source SI. The best fitting RS spectrum is 
overlaid and the residuals (Data— Model in counts) of the fit are shown in the lower panel. We 
marked the elements responsible for major line complexes. The inset displays the confidence 
regions in the Nh—T plane (68.3%, 90.0%, and 99.0% confidence levels; the best fit is marked 
by the plus). 




Fig. 7. — ACIS-S3 spectra of the point sources S2-S6. Overlaid are the best fits for the 
following models. S2 and S3: power law, S5 and S6: RS, S4: black body (see Table 2). 
The residuals of these fits (Data— Model in counts) are shown in the lower panels. The insets 
display N H — T confidence regions at 68.3%, 90.0%, and 99.0% levels; the plus indicates the 
best fit. Panels showing spectra and residuals are all on the same scale. 
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Fig. 8. — Azimuthally-averaged X-ray surface brightness profiles for different energy ranges, 
all on the same logarithmic scale. The curves are best exponential fits for radii < 700 pc. 
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Fig. 9. — Radial profile of the three dimensional volume density distribution of the hot gas. 
Best fits of the following functions are overlaid. Solid line: exponential, dotted line: power 
law, dashed line: power law with an index of /3 — 2. The inset of this figure is an example 
of how a profile using six shells is calculated (see Sect. 2.3 for details). 
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Segment No. 

Fig. 10. — Azimuthal dependence of the X-ray and Ha features in NGC3077. The upper 
panel shows the numbering of the individual segments displayed as black polygons. For 
illustration, the regions defined in Fig. 11 are overlaid in grey. Lower panel: The intensity 
per unit area of the total X-ray band and the Ha image within the corresponding segments 
in arbitrary units. 




Fig. 11. — Definition of the regions Rl to R7 overlaid on an Ha image. 
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Fig. 12. — a) The background-subtracted spectrum of the total diffuse emission of NGC 3077 
(Rl). Identified point sources were previously removed, b) Confidence regions for the two 
different plasma models (RS and MeKaL) in 68.3%, 90.0%, and 99.0% intervals. Solid: RS 
model (best fit visualized by a circle), Dotted: MeKaL plasma (plus). The corresponding 
best fitting spectra of the models are overlaid to the measured spectrum in panels c) and 
d), as well as the elements responsible for prominent line complexes. 
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Fig. 13.— Nh — T Confidence regions (68.3%, 90.0%, and 99.0%) of RS and MeKaL plasma 
models for the individual regions R2 to R7. The plus (circle) refers to the best fit of the 
model shown in solid (dotted) contours. All plots are on the same scale. First row: Both 
sets of contours are from RS plasma models. The dotted ones are for the overall diffuse 
emission of Rl, the solid ones for the individual bubbles. Second row: Like the first row, 
but for the MeKaL models. Third row: Overlay of the MeKaL models (dotted) on the 
RS models (solid) for each region (for the Rl, the diffuse emission of the entire galaxy, see 
Fig. 12). 
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Fig. 14. — Hardness ratio plot of the point sources in the Chandra ACIS-S3 field of view 
which were detected by both wavdetect and celldetect source detection algorithms. Open 
Squares: The point sources within the diffuse X-ray emission of NGC 3077. The sources 
are marked according to Fig. 1 and Table 2. Filled Circles: Background Sources. Objects 
which are identified as stars are plotted as filled triangles. Power law models are overlaid 
on the data. Solid lines: The photon indices 7 are kept fixed, while the absorbing column 
densities vary from to 10 22 cm~ 2 . Dotted lines: For a fixed absorbing column density Njj 
the photon indices change from to 4. Column densities are given in units of 10 21 cm -2 . The 
shaded region marks the positions of the regions Rl to R7. Thermal plasma models with 
temperatures below ~ 15 x 10 6 K occupy the region in the upper left corner - well above the 
power law models. 
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Fig. 15. — Total X-ray spectrum of NGC 3077. Overlaid are the best fits to the point sources 
and to the diffuse emission regions, the latter modeled by a): RS models and b): MeKaL 
models. Rl corresponds to the total diffuse X-ray emission of NGC 3077. 
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Fig. 16. — LogN(> S)-LogS plots of point sources in the G01 soft and hard bands within 
the ACIS-S3 field of view (binned to 1" pixel size). We excluded sources found within the 
hot outflow of NGC3077 as well as optically identified stars (corresponding to objects 3, 4, 
and 7 in Fig. 1 and Table 1). The data are compared to the results from G01 for the CDFS 
(solid lines), where the shaded regions mark the error of their fit. The limiting fluxes 
(S/N < 3) are indicated by dashed lines. 
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Fig. 17. — H I column density profile along the north-south direction of NGC3077 through 
the H-band peak (the center of NGC3077). For a comparison, we also show a profile of the 
diffuse X-ray emission along the same direction. The latter has been scaled to match the 
H I peak intensity. 



Table 1. Source list of the ACIS-S3 point sources (exposure time: 53.4 ks). 



No. 


Algorithm RA (J2000) a 


DEC (J2000) a 


Total Counts 


Soft Counts 


Medium Counts Hard Counts 


HR1 


HR2 


(1) 


(2) (3) 


(4) 


(5) 


(6) 


(7) (8) 


(9) 


(10) 



*1 


w 


10 


02 


45.2 


68 


46 


18.9 


11.0 ± 3.5 


0.0 ± 0.0 


0.0 ± 0.0 


11.4 ± 3.5 


— 1.00 ± 0.00 


— 1.00 ± 0.00 


2 




10 


02 


54.2 


68 


44 


57.5 


68 5 + 8 4 


6.7 ± 2.6 




50 9 + 7 2 


— 80 + 09 


— n 49 + n 11 


*3 b 


w /c 


10 


02 


56.0 


68 


47 


06.9 


333.8 ± 18.4 


115.1 ± 10.8 


191.8 ± 13.9 


26.9 ± 5.4 


—0.31 ± 0.06 


+0.84 ± 0.03 


*4c 


w/c 


10 


02 


56.4 


68 


47 


09.5 


1273.6 ± 35.7 


314.3 ± 17.7 


823.7 ± 28.7 


135.5 ± 11.8 


—0.51 ± 0.03 


+0.79 ± 0.02 


5 


w /c 


10 


02 


56.7 


68 


44 


30.5 


140.5 ± 11.9 


6.6 ± 2.7 


39.7 ± 6.3 


94.2 ± 9.7 


—0.91 ± 0.06 


—0.34 ± 0.08 


*6 


w 


10 


03 


00.0 


68 


44 


07.8 


5.6 ± 2.5 


0.0 ± 0.0 


0.0 ± 0.0 


5.7 ± 2.5 


— 1.00 ± 0.00 


— 1.00 ± 0.00 


*7 




10 


03 


08.2 


68 


41 


05.8 


116 4 + 10 9 


39 4 + 6 3 


40.9 ± 6.4 


36.1 ± 6.1 


—0 32 + 1 1 


+0 38 + 09 

i yj * '-J o _i_ yj . w u 


8 


w /c 


10 


03 


08.4 


68 


44 


08.6 


28.2 ± 5.4 


0.8 ± 1.0 


9.9 ± 3.2 


17.5 ± 4.2 


—0.94 ± 0.13 


—0.24 ± 0.19 


9 


w/c 


10 


03 


09.9 


68 


43 


13.3 


14.2 ± 3.9 


0.0 ± 0.0 


0.0 ± 0.0 


14.5 ± 3.9 


— 1.00 ± 0.00 


— 1.00 ± 0.00 


10 d 


w /c 


10 


03 


10.5 


68 


41 


50.2 


43.2 ± 6.6 


4.8 ± 2.2 


10.0 ± 3.2 


28.4 ± 5.4 


—0.78 ± 0.13 


—0.31 ± 0.15 


11 


w /c 


10 


03 


12.3 


68 


43 


18.9 


10.5 ± 3.3 


0.0 ± 0.0 


0.9 ± 1.0 


9.8 ± 3.2 


— 1.00 ± 0.05 


—0.83 ± 0.18 


12 


w /c 


10 


03 


20.6 


68 


45 


07.1 


19.3 ± 4.5 


3.7 ± 2.0 


7.0 ± 2.6 


8.7 ± 3.0 


—0.62 ± 0.24 


+0.10 ± 0.23 


*13 


w /c 


10 


03 


20.8 


68 


41 


40.0 


221.0 ± 14.9 


30.8 ± 5.6 


71.9 ± 8.5 


118.3 ± 10.9 


—0.72 ± 0.06 


—0.07 ± 0.07 


*14 


w 


10 


03 


21.4 


68 


46 


33.4 


6.5 ± 2.6 


1.8 ± 1.4 


2.0 ± 1.4 


2.7 ± 1.7 


—0.45 ± 0.45 


+0.17 ± 0.40 


15 


w /c 


10 


03 


22.0 


68 


45 


03.2 


12.1 ± 3.6 


0.8 ± 1.0 


4.9 ± 2.2 


6.4 ± 2.6 


—0.87 ± 0.26 


—0.06 ± 0.29 


*16 


w /c 


10 


03 


23.0 


68 


40 


31.7 


32.4 ± 5.8 


3.5 ± 2.0 


3.9 ± 2.0 


25.0 ±5.1 


—0.78 ± 0.13 


—0.54 ± 0.15 


*17 


w /x 


10 


03 


26.6 


68 


39 


22.5 


9.4 ± 3.3 


0.0 ± 0.0 


1.8 ± 1.4 


7.9 ± 3.0 


— 1.00 ± 0.11 


—0.63 ± 0.26 


*18 


w 


10 


03 


26.8 


68 


41 


29.6 


3.8 ± 2.0 


0.0 ± 0.0 


0.0 ± 0.0 


3.9 ± 2.0 


— 1.00 ± 0.00 


— 1.00 ± 0.00 


*19 


w 


10 


03 


27.2 


68 


41 


59.3 


5.8 ± 2.4 


0.0 ± 0.0 


1.0 ± 1.0 


4.9 ± 2.2 


— 1.00 ± 0.13 


—0.66 ± 0.31 


*20 


w /c 


10 


03 


29.0 


68 


44 


25.3 


11.5 ± 3.5 


0.0 ± 0.0 


1.9 ± 1.4 


9.7 ± 3.2 


— 1.00 ± 0.09 


—0.67 ± 0.22 


*21 


w /c 


10 


03 


29.9 


68 


44 


22.8 


18.5 ± 4.4 


4.8 ± 2.2 


3.9 ± 2.0 


9.7 ± 3.2 


—0.48 ± 0.25 


—0.05 ± 0.24 


22 


w/c 


10 


03 


30.0 


68 


40 


41.8 


32.5 ±5.8 


4.7 ±2.2 


3.9 ±2.0 


24.0 ±5.0 


-0.71 ±0.14 


-0.47 ±0.16 


*23 


w/c 


10 


03 


31.0 


68 


41 


54.2 


20.3 ±4.6 


1.8 ±1.4 


1.0 ±1.0 


17.6 ±4.2 


-0.82 ±0.14 


-0.73 ±0.16 


*24 


w/c 


10 


03 


32.9 


68 


46 


17.0 


31.0 ±5.7 


7.8 ±2.8 


7.9 ±2.8 


15.3 ±4.0 


-0.50 ±0.19 


+0.01 + 0.18 


25 


w 


10 


03 


33.4 


68 


40 


05.2 


3.9 ±2.0 


0.0 ±0.0 


1.0 ±1.0 


2.9 ±1.7 


-1.00 ±0.22 


-0.49 ± 0.44 


*26 


w/c 


10 


03 


41.8 


68 


41 


29.7 


150.9 ± 12.3 


15.7 ±4.0 


32.8 ±5.7 


102.4 ± 10.1 


-0.79 ±0.06 


-0.36 ±0.08 


*27 


w 


10 


03 


43.7 


68 


44 


33.5 


11.1 ±3.5 


0.8 ±1.0 


0.0 ±0.0 


10.4 ±3.3 


-0.86 ±0.17 


-0.86 ±0.17 


28 


w/c 


10 


03 


45.9 


68 


41 


55.7 


21.8 ±4.8 


1.7 ±1.4 


1.9 ±1.4 


18.3 ±4.4 


-0.84 ±0.14 


-0.67 ±0.16 


29 


w/c 


10 


03 


49.0 


68 


42 


56.3 


19.5 ±4.5 


0.9 ±1.0 


3.9 ±2.0 


14.7 ±3.9 


-0.91 ±0.14 


-0.51 ±0.20 


*30 


w/c 


10 


03 


49.1 


68 


45 


51.1 


29.7 ±5.6 


1.8 ±1.4 


4.8 ±2.2 


23.1 ±4.9 


-0.88 ±0.12 


-0.56 ±0.16 


*31 


w 


10 


03 


53.0 


68 


42 


22.9 


13.7±3.7 


0.0 ±0.0 


3.0 ±1.7 


10.8 ±3.3 


-1.00 ±0.10 


-0.57 ±0.22 


32 


w 


10 


03 


54.6 


68 


41 


24.8 


15.9 ±4.5 


3.3 ±2.0 


3.3 ±2.0 


9.4 ±3.5 


-0.59 ±0.27 


-0.18 ±0.28 


33 


w 


10 


03 


55.1 


68 


41 


08.6 


6.4 ±2.6 


0.8 ±1.0 


0.0 ±0.0 


5.6 ±2.5 


-0.75 ±0.29 


-0.75 ±0.29 


*34 


w/x 


10 


03 


57.2 


68 


40 


44.4 


11.2 ±3.8 


0.0 ±0.0 


2.6 ±1.7 


9.9 ±3.3 


-1.00 ±0.11 


-0.58 ±0.24 


*35 


w/c 


10 


03 


58.5 


68 


43 


41.6 


265.2 ± 16.4 


37.3 ±6.2 


82.6 ±9.1 


145.4 ± 12.1 


-0.72 ±0.06 


-0.10 ±0.06 


*36 


w/c 


10 


04 


02.8 


68 


42 


19.4 


17.8 ±4.7 


0.0 ±0.0 


1.5 ±1.4 


17.3 ±4.5 


-1.00 ±0.04 


-0.84 ±0.14 



Table 1 — Continued 



No. 


Algorithm 


RA (J2000) a 


DEC (J2000) a 


Total Counts 


Soft Counts 


Medium Counts 


Hard Counts 


HR1 


HR2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


*37 


w/x 


10 04 16.2 


68 44 35.4 


25.6 ±5.4 


2.2 ± 1.7 


3.7±2.0 


19.7 ±4.7 


-0.83 ± 0.15 


-0.54 ±0.18 


*38 


w/x 


10 04 29.4 


68 43 36.6 


60.6 ± 8.8 


9.0 ±3.7 


15.3 ±4.2 


36.3 ± 6.8 


-0.70 ± 0.14 


-0.20 ± 0.14 



a Right ascension is given in h m s, Declination in ° ' " 

b Foreground star ADS 7611 A (binary system with Object No. 4, HD 86677) 
foreground star HD 86677 (BD+69 552) 

d 8" distance to star J100311.93±684148.3 (Sakai & Madorc 2001) 



Table 2. Spectral fitting results of the point sources in NGC 3077. 



Parameter 


Unit 


SI 


S2 


S3 


S4 


S5 


S6 


RA (J2000) 




10 h 03 m 18 ! ?8 


10 fc 03 m 19 ! rl 


l0 h 0Z m WA 


10 ,, 03 m 17 ! :8 


10 h 03 m 17?9 


10 ,l 03 m 18?3 


DEC (J2000) 




68°43'56.'4 


68°44'01'.'4 


68°44'02'/3 


68°44'16'/0 


68° 43' 57'.' 3 


68°44'03.'8 


Photons 


[cts] 


133 ± 12 


114 ±11 


119 ± 11 


37 ±7 


17 ±4 


17 ±4 


Count Rate 


[lCT 3 cts s- 1 ] 


2.50 ±0.21 


2.13 ± 0.21 


2.23 ± 0.21 


0.69 ± 0.13 


0.32 ± 0.07 


0.32 ±0.07 


RS Model 


N H 
T 

Norm. a 

rpabs 
r X 

Fx 
Lx 


[10 21 cm- 2 ] 
[10 6 K] 

[10- 15 crg cm" 2 s- 1 ] 
[10 _15 erg cm -2 s -1 ] 
[10 37 crg s- 1 ] 


9-50lJ'g2 
9 10+ ' 87 

y - lu -0.85 
o oo + l.OO 
J - ZO -0.76 

10.48±I;i? 
96.58™ 
14.98ll;g| 


13.67 +3 - 7 7 
102 +i59 

luz -41.2 
n oc+0.82 
z - oo -0.19 

39.88+. ^-a? 
63.16t^ 25 


pr oc + 1-50 

O-^O-l.lO 
Q 9 OO+84.90 

2 50+ ' 47 
Z.OU_ 23 

39-55^ 2 6 g 3 
54.89li 1 7 I 3 

8-5i±i;!^ 


2-78t?; 5 28 
65+ ' 22 
o 04 +64 
ou - r 207 
q oo+5'22 

Z -*J_1.97 


7.86+°'^ 

06+ ' 11 

u ' UD -0.02 

1 36+ 4 ' 46 
1 "- >o -0.92 


R n1 +2.42 

6.01I 204 

*-" 3 ' -1.42 

0.58+ 4 ;^ 

1 7 09 +181 

2 65+ 28 ' 
Z - D0 -1.45 



Power Law Model 





[10 21 cm- 2 ] 5 


5 . 55 


14.48+1.3 7 


q qq+1-22 
O ' OO -1.00 


n r^+0.99 
u.oo_q 45 


4.08+ 3 ' 7 3 


1 QQ+3.15 
±u -°°-2.30 


7 




o 9^+0.43 


i-65±8:ti 


1 00+ ' 26 
l.UU_ 25 


5 03+ 1 ' 44 


4 p;4+ 2:i3 
*- 04 -l.l)p 


o 97+1.73 

-0.60 


Amplitude 


[IO" 6 ] 


15.57±2 ;i | 


9.09+^ 


q 7 7+l.ll 

°. ' ' -0.79 


S2+ 013 


9 q^+5.85 


19.08± 38 85 4 


rpabs 
1 X 


[10- 15 crg cm - 2 s^ 1 ] 




41 04+ 129 


oo.ua_28.67 


q cq+76.43 


i kq+27.41 
32 


1.41™ 


Fx 


[10 _15 erg cm - 2 s^ 1 ] 


124+ 51 9 7 


71 5fi+ 161 
<1.0D_ 31 45 


65.14ti;? 3 


14 27+ 101 
i4 '^' -9.58 


62.39l 2 « 5 1 ° 5 


49000+ 116000 
4yuuu -37000 


Lx 


[10 37 crg s- 1 ] 


19.25l^ 1 


11 10+ 24 ' 95 
11.1U_ 4 gg 


10.10^; 5 5 2 


9 9-1 +15.72 
z - zi -1.49 


9-67l 7 1 6 1 2 


7^70+180000 



Black Body Model 





[10 21 cm- 2 ] 


0.98t°;f 


7.26+J; 7 ? 


0.581°;°° 


„ qO+0.61 

U -38_ .38 




4.36tli 3 


T 


[10 6 K] 


5.88+° 4 l 


12.04+ / 3 ; 7 


13.62±iii 


94+° ll 
u - y4 -0.12 


q r-r-+0.29 
°. oo -1.00 


i-8i!8:!S 


Amplitude 11 




2 19+ 1 ' 24 
z - ly -0.65 


29+ ' 23 
u ' zy -0.12 


20+ 09 


9n9 9+6020 
^ u ^-1820 


i OO + 31.08 
l-88_ 01 


234 j; 287 


T^abs 

b x 


[10- 15 erg cm -2 s-i] 


12.95t 7 7 5 f 


,o fi1 +60.98 

q/i nc+65.87 
o^.zo_22.43 


3^.49_ 20 .4 9 


9 „o + 3(i.H2 
A °°-2.83 


i fi o+37.34 
i - D ' :5 -1.42 


i qo + 14.20 
1.38-1.20 


Fx 


[10- 15 erg cm -2 s -1 ] 


14.92l 7 6 8 f 


, 9 10+40.30 
38.18_20.29 


t- 09+33.88 
°- yZ -5.63 


1.69+ 38 2 78 


13.04l 47 5 21 


Lx 


[10 37 crg s- 1 ] 


9 01 +2.59 


t- qi +10.22 
J - 01 -3.48 


5-92±|;f| 


n O9 + 05-25 
u - yz -0.87 


26+ 6 ' 01 
U - ZD -0.19 


2.02t 7 ; 32 



a The normalization is given in units of 10- 5 K, where K is IO-1 4 (4ttD 2 y 1 J n e n p dV; Da is the angular size distance to the source 
[cm], n e and n p are the electron and proton densities [cm" 3 ]. 

b The amplitude is given in units of IO -5 A, where A = 2ivc^ 2 h^' A {R/d) 2 = 9.884 X 10 31 (R/ti) 2 ; c is the speed of light [cm s _1 ], h 
is Planck's constant [keVs], R is the radius, and d the distance to the source. 



Table 3. Best fit parameters of spectral RS and MeKaL models to the regions Rl to R7. 



Parameter 


Unit 


Rl (Total) 


R2 


R3 


R4 


R5 


R6 


R7 


Area 


[arcsec 2 ] 


8463 


268 


256 


294 


341 


228 


220 


Net Counts 


[cts] 


729 ± 46 


142 ± 14 


104 ± 12 


43 ±9 


58 ± 11 


71 ± 10 


58 ± 10 


Net Count Rate 


[KT 3 cts s- 1 ] 


13.7 ±0.4 


2.7 ±0.0 


1.9 ±0.2 


0.8 ±0.2 


1.1 ± 0.2 


1.3 ±0.2 


1.1 ±0.2 


RS Best Fit Parameters 


N„ 
T 

Norm. a 

J?abs 
r X 

Fx 
Lx 


[10 21 cm- 2 ] 
[10 6 K] 

[10- 15 crg cm- 2 s- 1 ] 
[10- 15 crg cm- 2 s" 1 ] 
[10 37 crg s- 1 ] 


D - O1 -0.32 
2 03+ 01 

31.94±||-I 

2900+™° 
450+mo 


6 54+1-00 

2 02 +a4() 
z - u -0.28 
2Q 41+43.42 

fi 59+55.7 
D - oz -5.93 
ooq + 1140 
663 -5 13 
1fl o + l77 


4 14+2.00 

9 7R +0.80 
z - ,D -0.67 

1 76+ 11 - 76 

1 -'°-1.04 

Q 0O+81.7 

3 11 

45.14+1* 
7.00± 49 4 6 5 


7 R7 +10.33 
'• D '-1.99 
1 9 7 +1.10 
1 - z '-0.37 

19411 22 

l-98li 8 9 9 7 
2300+ 5i7 ° 

356+ 801 
*J°_306 


9.77+^ 

2-68l 4 2 5 6 8 7 
194+. 47 f 
30.05+J 4 ;; 4 


5 07+1- 56 
°- u '-1.34 

2.3l±°, : f 7 

4-06™ 
3.08± 49 3 7 

98.00+ 4 q\ 

15.20+™ 


9 41+0.86 

M -*l-1.01 

1 59+°' 25 
1- OM -0.16 
+ 169 
10, -136 

2.68+1% 

9700+3960 
z,uu -2390 
41 o+6114 
* io -371 



McKaL Best Fit Parameters 



N„ 


[10 21 cm- 2 ] 


5 01+°' 72 
°- ui -0.41 


4-61±?Il 


4.68+ij 3 


7.58t°;? 7 




0-4to'.4 3 


9 99+0.59 

y - zy -0.99 


T 


[10 6 K] 


9 04+0. 08 


2-601°;!° 


2-48t°;S 


1 9C+0.23 

l- zo -0.11 


9 1O+0.45 
^•1»_0.37 


4.88t°;f 2 


-, OO+0.29 
1- JO -0.03 


Norm. a 




46.07±||;? 


5 56+ 6 ' 24 
d - JU -6.32 


/17+22.57 

<i -*' -1.65 


2=0+81.5 

80.2 


2.90+f 5 °; 6 


15+°- 25 
u - lo -0.03 


ooo+79.8 
oou_ 78 7 


Triads 

b x 


[10- 15 crg cm" 2 s-i] 


32.82+ 4 ^ 


6-58l 45 i 


S 4fi+ 82,7 


2-00+" 3 g 


2 74+11 5 
^•'+-2.49 


9^+7.98 
°- ZJ -2.03 


2.70t 9 / 7 j l 


Fx 


[10- 15 crg cm- 2 s-i] 


o 7 o+i530 
y ' u -69.4 


124+J 9 2 


75.96± 4 28 2 


2510™ 


58.24t 725 5 


4-i6t 7 ; 1 13 


3940li|?° 


Lx 


[10 37 crg s-i] 


1 5n+ 238 

iou -10.8 


19.25l 3 2:4 


11.78±|^| 


389+J 42 


9-03±^ 5 2 o 


0-65tJ;i° 


fil i+206 
Dli -296 



a See footnote a in Table 2 



Table 4. Derived parameters of the hot gas from the best RS and MeKaL fits. 



Parameter 


Unit 


Rl (Total) 


R2 


R3 


R4 


R5 


R6 


R7 






Derived Parameters of the Hot Gas a 


(RS/McKaL) 








deq 


[PC] 


1811 


322 


315 


337 


306 


297 


292 


mean line of sight 


[PC] 


1207 


215 


210 


225 


204 


198 


195 




[cm- 3 ] 


0.05/0.03 


0.30/0.13 


0.08/0.11 


0.71/0.82 


0.19/0.10 


0.13/0.02 


0.80/1.16 


EMlxiUv)- 1 ] 


[cm -6 pc] 


3.0/1.1 


19.4/3.6 


1.3/2.5 


113.4/151.3 


7.4/2.04 


3.4/0.1 


124.8/262.4 




[10 B Kcm- 3 ] 


2.0/1.4 


12.1/6.8 


4.4/5.5 


18.0/20.5 


6.8/4.4 


6.0/2.0 


25.4/32.0 




[10 4 M Q ] 


384.2/230.5 


13.0/5.6 


3.2/4.4 


35.2/40.7 


7.0/3.7 


4.4/0.7 


25.8/37.4 




[10 52 erg] 


384.1/265.7 


12.9/7.2 


4.4/5.4 


22.0/25.0 


6.2/4.0 


5.0/1.6 


20.2/25.4 


ui[xr°- 5 /„ a5 ] 


[Myr] 


27.0/55.9 


4.0/11.8 


19.9/14.5 


2.0/2.0 


6.5/14.0 


10.4/77.9 


1.5/1.3 


M coo i 


[Meyr- 1 ] 


0.142/0.041 


0.033/0.005 


0.002/0.003 


0.176/0.204 


0.011/0.003 


0.004/0.0001 


0.172/0.288 


< ""hot > 


[kms^ 1 ] 


300/380 


320/360 


370/350 


250/250 


300/330 


340/480 


280/260 



Predicted Parameters' 3 



Notation of Shells 


in M98 


D 


A 


G 


B 


J 




[cm -3 ] 


0.45 


0.60 


1.09 


1.55 


2.34 


^mech 


[lO^ergs- 1 ] 


2.0 


2.6 


3.9 


1.6 


25.7 


Lx 


[lO^crgs- 1 ] 


11.3 


38.1 


82.4 


40.0 


667.8 


..pro j 
u exp 

^shell 


[km s" 1 ] 
[Myr] 


106 
1.7 


55 
8.6 


51 
10 


40 

9.4 


67 
8.2 



a See Sect. 2.4. 
b See Sect. 3.2.1. 
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Table 5. Best fitting MeKaL models for the probable supernova remnants. 



Parameter 


Unit 


SI 


S5 


S6 


MeKaL Model 


N H 
T 

Norm. a 

j?abs 

r X 
Fx 
L x 


[10 21 cm" 2 ] 
[10 6 K] 

[10 -15 erg cm -2 
[10 -15 erg cm -2 
[10 37 erg s- 1 ] 


9-iv±i;^ 

q 07+ 113 

q OO + 1-05 
°- zo -0.85 
s-l] 10.44t 8 4 .g 

s-i] 85.75+^.31 
13-30±t;gg 


1 18+ 1 ' 94 

7.36± 1 1 °I 

o.o8±g-;s? 

i 07+O6.35 
1 - 3 '-0.88 
9 9s +7.20 
^•^°-0.46 

U - JO -0.07 


u - z '-0.10 
1 40+ 20 - 84 
7 c 9 +52.68 
1 - J -3.30 
1 17 +8.17 
J-- 1 ' -0.51 



a See footnote a in Table 2 



Table 6. Comparison^ NGC 3077, M 81, and M 82. 



Parameter 


Unit 


NGC 3077 a 


M81 


M81 b 


M81 b ' c 


M 82 M 82 b 


Reference 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) (8) 


(9) 


Total L<£ tai 


[10 38 erg s" 1 ] 


49/19 


283 




97 


351 196 


1,2,3 


. . r point 

Point sources L x 


[10 38 erg s" 1 ] 


4 (8%/21%/27%) 


210 (74%) 


27 (28%) 


88 (91%) 


165 (47%) 10 (5%) 


1,2,3 


Diflfusc emission L% f fuse 


[10 38 erg s" 1 ] 


45/15 (92%/79%/63%) 


70 (25%) 


70 (72%) 


9 (9%) 


186 (53%) 186 (95%) 


1,2,3 


Total Mass Mtot 


[10 1,J M ] 


1.2 




12 




3.1 


4,4,4 


Hi Mass M HI 


[10 8 M ] 


1.3 




27 




8.8 


5,6,6 


Blue Luminosity d LB 


[10 s L Q ] 


14 




227 




59 


7,7,7 


Infrared Luminosity Lju 


[1O 8 L ] 


3.2 




2.5 




207 


8,8,8 


Ho Luminosity Ljj a 


[10 38 erg s" 1 ] 


81 




540 




4962 


5,9,10 


SFR C 


[Meyear- 1 ] 


0.06 




0.4 




4 





Normalized Luminosities 



L tot a l/ MHI 


[10 30 erg s"VM ] 


37.7/14.6 


10.5 


3.6 


39.9 




22.3 


L p ° int /M HI 


[10 30 erg s-VM ] 


3.1 


7.8 


1.0 3.3 


18.7 




1.1 


L™ ffuse /M HI 


[10 30 erg s-'/Mq] 


34.6/11.5 


2.6 


2.6 0.3 


21.1 




21.1 


L B /M HI 


[Lq/Mq] 


10.8 




8.4 




6.7 




Lir/M hi 


[Lq/Mq] 


2.5 




0.09 




23.5 




L Ha /M HI 


[10 30 erg s- 1 /Mq] 


62.3 




20.0 




564 




SFR/M HI 


[lO-Sycar- 1 ] 


0.05 




0.01 




0.45 




M hotga s/M HI 




0.030/0.018 








0.15 




E th /M HI 


[10 46 erg/M ] 


3.0/2.0 








41 




M coo i/M HI 


[ 10 -10 yr -l] 


10.9/3.1 








2.3 





1 For all three galaxies we adopt a common distance of 3.6 Mpc. All numbers taken from the literature are corrected respectively. 
a for both models: RS/MeKaL; the relative luminosities given in parenthesis refer to RS/MeKaL/Count Rate 
b without the strong point— like nuclear source 
c corrected for unresolved point sources 

d using the extinction— corrected magnitude rrfig^ and a solar blue magnitude of 5.50. 

e Calculated from the Ho luminosities via SFR = L Ha /(1.26 X 10 41 erg s" 1 ) M yr" 1 (Kennicutt, Tamblyn, & Congdon 1994). 

References. — (1) this work; (2) Immler & Wang 2001; (3) S97; (4) estimated by assuming an M/L » 9 by Brouillet et al. 1991; (5) W02; (6) Appleton, 
Davies, & Stephenson 1981; (7) dc Vaucoulcurs ct al. 1991; (8) Yun 1999; (9) Grccnawalt, Walterbos, Thilker, & Hoopcs 1998; (10) Young, Klcinmann, & 
Allen 1988. 



